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I Abstract 

Relativistic domain walls are studied in the framework of a polynomial ap- 
proximation to the field interpolating between different vacua and forming 
the domain wall. In this approach we can calculate evolution of a core and of 
r-| ! a width of the domain wall. In the single, cubic polynomial approximation 

used in this paper, the core obeys Nambu-Goto equation for a relativistic 
membrane. The width of the domain wall obeys a nonlinear equation which 
is solved perturbatively. There are two types of corrections to the constant 
zeroth order width: the ones oscillating in time, and the corrections directly 
related to curvature of the core. We find that curving a static domain wall 
is associated with an increase of its width. As an example, evolution of a 
toroidal domain wall is investigated. 
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1 Introduction 



Recently one observes rapidly growing interest in time evolution of topolog- 
ical defects in 3+1 dimensional space-time. This subject is important for 
many branches of physics. Without attempting to present here a complete 
list, let us mention vortices in superconductors |jl|] and in superfluids 0, de- 
fects in liquid crystals 0, magnetic domain walls 0], cosmic strings 0, [|], 
and a flux tube in QCD 0. While evolution of topological defects in 1+1 
dimensional space-time has been rather well understood, in 3+1 dimensions 
relatively little is known, and the problem is actually a formidable one. 

The type of equations from which one attempts to calculate evolution of 
topological defects depends in an essential manner on the physical context. 
In condensed matter physics one uses e.g. diffusion type equations |^ or 
nonlinear Schrodinger type equations 0. For cosmic strings in a negligi- 
ble gravitational field, or for the particle physics flux-tubes, one should use 
Poincare invariant wave equations. 

Our paper is devoted to dynamics of domain walls governed by a Poincare 
invariant wave equation. In this case, several analytical approaches have 



been made, see e.g.H,[ttn],|ll|,[12] as well as numerical calculations, see, e.g. 



[|r^,[n|. They have given insights into the dynamics of domain walls, and 
have revealed the richness and intricacy of it. Our motivation for studying 
the dynamics of relativistic domain walls is of rather mathematical character 
- relative simplicity of the pertinent field equations makes them a convenient 
testing ground for new methods of calculating the evolution of topological de- 
fects. Nevertheless, direct physical applications are also possible. Relativistic 
domain walls appear in a field-theoretic approach to cosmology, and in par- 
ticle physics (e.g. the surface of a quark bag can be regarded as a domain 
wall). Moreover, the polynomial approximation we develop in this paper can 
also be applied to non-relativistic domain walls observed in condensed matter 
physics - no significant changes are required. 

Till now, the main line of the analytical approaches to the description of 
dynamics of a domain wall has been to reduce the initial, 3+1 dimensional 
field-theoretical system to an effective theory of a classical, relativistic mem- 
brane. This approach, called the effective action method, is very appealing 
conceptually, and it is correct in principle. However, because of the complex- 
ity of the pertinent field-theoretic equations, it is rather difficult to carry out 
the necessary calculations without shortcuts, which in turn introduce some 
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uncertainty about the final result. For a critical discussion of the effective 



action method we refer the reader to the paper |T^. For the most recent 
application of the effective action method to domain walls, see [|r^. 

Our opinion is that at the present stage of the subject one should develop 
and refine various methods of investigations of the dynamics of the domain 
walls. It seems that because of enormous complexity of the full, 3+1 dimen- 
sional, nonlinear field-theoretic dynamics involved, it is long time until we 
can analytically or numerically calculate evolution of a generic domain wall 
without any difficulty. 

In the present paper we generalize the method of analysis of evolution of 
relativistic domain walls proposed in . The main characteristic feature of 
this approach is a simple, approximate polynomial Ansatz for the field inside 
the domain wall, while outside of it the field has the exact vacuum values. 
The coefficients of this polynomial are calculated from the field equation and 
from boundary conditions. Actually, this approximation can be regarded 



as an application of splines, |18|. In the paper this approach has been 
applied only to cylindrical and spherical domain walls. In the present paper 
we apply the polynomial approximation to generic smooth domain walls - 
this is presented in Section 2. We approximate the field inside the domain 
wall by a cubic polynomial in a transverse, co-moving coordinate. We obtain 
Nambu-Goto type equation for a core of the domain wall, and a nonlinear 
equation for the width of the domain wall. Also, a perturbative scheme for 
solving the latter equation is presented. It is based on dividing the width 
into two components: the one oscillating with characteristic frequency given 
by the mass of the scalar field, and the other one directly related to curvature 
of the core. 

In the next Section we rewrite the Nambu-Goto equation, and a formula 
for the width of the domain wall, in terms of local curvature radia and velocity 
perpendicular to the core. This gives a rather nice and useful insight into the 
dynamics of the domain wall. We show that a force acting on a small piece 
of the core can be regarded as being due to a surface tension. We also show 
that curving a static domain wall is associated with an increase of its width. 

In Section 4 we calculate the energy of the domain wall. We find that the 
energy density depends on the curvature of the core, and that nonuniformities 
of the width increase the energy density. We also find, rather surprisingly, 
that curving the domain wall seems to decrease its energy. These results, as 
well as all others in this paper, are obtained for slightly curved domain walls; 
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our approximations break down when the curvature radia become comparable 
with the width of the domain wall. 

Next, in Section 5, we consider as an example evolution of a toroidal 
domain wall. We find two types of evolution of this domain wall. 

Section 6 contains a suggestion how to improve our approximate Ansatz 
for the scalar field, and also some ending remarks. 

In the Appendix we discuss accuracy of the polynomial approximation. 



General formalism for evolution of the do- 
main wall in the polynomial approximation 



In this subsection we would like to generalize the formalism presented in IT? 
where only cylindrical and spherical domain walls were considered, and to 
present the corresponding, approximate solution of the field equation. 

We will investigate domain walls in the well-known model, see e.g. [^], 
involving only a single, real scalar field $ with the following Lagrangian 

L = -lv,.d'^<^d^^-^{<^'-^r, (1) 

where {r]^^) = diag{—l, 1, 1, 1), and A, M are positive constants. The corre- 
sponding field equation is 

- 2A($2 -)$ = 0. (2) 

4A 

The vacuum values of the field $ are equal to ±$o, where $o = M/2\^. 
The domain wall arises if at a given time the field is equal to one of the 
two vacuum values in some region of the space, and is equal to the other 
vacuum value in the complementary part of the space, except for the border 
layer between the two regions (the domain wall), where the field smoothly 
interpolates between the vacuum values. It is clear that at each instant of 
time the field $ vanishes somewhere inside the border layer. We assume 
that the locus of these zeros is a smooth surface S. We shall call it the core 
of the domain wall. The well-known example of the domain wall, with the 
static core given by the (x^, x^) plane, is given by the following exact, static 
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solution of Eq.(2) 



$ = $0 tanh{—), (3) 

where Iq = M~^. The width of this domain wall is of the order Iq, and energy 
density is exponentially localised around the {x^,x^) plane. 

For a generic domain wall, space-time parametrisation of the world- volume 
S of the core (a 3-dimensional manifold embedded in Minkowski space-time, 
whose time slices coincide with S) can be chosen as follows 



( - \ 

\X%T,a\a') J 



(4) 



where r coincides with the laboratory frame time x^, and cr^, parametrise 
the core 5* at each instant of time. 

As usual, we introduce a special coordinate system (r, cr^, cr^, ^) in a vicin- 
ity of the world- volume S, co-moving with the core, 0. The new coordinates 
(r, cr^, cr^, ^) are defined by the following formula 

x^" = X^{T,(x\a^) +^n^'{T,a\a^), (5) 

where x'^ are Cartesian lab-frame coordinates in Minkowski space-time, and 
{n'^) is a normalised space-like four-vector, orthogonal to the S (in the co- 
variant sense) i.e. 

n^X^, = 0, = 1, (6) 

where a = corresponds to r; a = l,a = 2 correspond to o"^,(T^, and 
X'^ = dX^/dr, etc. The four- vectors X^^, X^-i, X0.2 are tangent to S. For 
points lying on the core ^ = 0, and the parameter r coincides with the lab- 
frame time For ,^ 7^ r is not equal to the lab-frame time x". The 
advantage of using the co-moving coordinates is that the world-volume S 
is described by the simple condition ^ = 0. Notice that the definition (5) 
implies that ^ is a Lorentz scalar. 

The next step is to write Eq.(2) in the new coordinates. It is convenient 
to introduce extrinsic curvature coefficients Kat and induced metrics gab on 
E: 

Kat = n,X%, gab = X^^„X^,,, (7) 
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where a,b = 0, 1, 2. The covariant metric tensor in the new coordinates can 
be readily calculated, and it can be written in the following form 



Gab 

1 



(8) 



where a, (3 — 0,1,2,3; a — 3 corresponds to the ^ coordinate; and 

Gab = Maeg"'Md,, Mac = 9ac-^Kac. (9) 

Thus, G^^ = 1, G^a = (a=0,l,2, as in (6)). It follows from formula (9) that 
\/—G, where G = det[Gai3], is given by the formula 

= detM, (10) 

where as usual g = det[gab]- detM can be explicitely evaluated in terms of 
Kab and gab'- 

detM = hir, a\ a\ g{T, a\ a^), (11) 

where 

KT,a\a\i) = l-^K:+^e{K:Kl-KK^)-y'eabce''fK2K',Kj. (12) 

In formula (11) we have noted that g can depend on r, a^,a^. Thus, 

V=G = ^ h. (13) 

For raising and lowering the latin indices of the extrinsic curvature coefficients 
we use the induced metric tensors gab, g"'^- The inverse metric tensor G'*^ 
is given by 

- Qab Q - 

1 



[G 



a/3l 



where 



G"^ = {M-^yg^diM-^)'^. 
Simple algebraic calculation gives exphcit formula for (M"^)"*^: 

{M-'r = I {^"11 - CK', + le{K',K', - i^W)] 



(14) 
(15) 



(16) 
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(this is just the matrix inverse to [Mab]; by definition it has upper indices). 
In general, the coordinates (r, cr^, cr^, ^) are defined locally, in a vicinity of the 
world- volume S. The allowed range of the ^ coordinate can be determined 
from the condition h > 0. Detailed discussion of the region of validity of the 



co-moving coordinates has been given in |T^. In the co- moving coordinates 



the field equation (2) has the following form 

1 ^-(v^ + hih 5,$) - 2A($^ - = 0. (17) 



The basic feature of our approach is the approximate cubic polynomial 
Ansatz for the $ field inside the domain wall. We assume that 

f +$0 for e>eo, 

$(r,a^a^O= A^ + ^Be + ^^Ce for -6<e<eo, (18) 
[ -$o for e < -6- 

Here E,OiE,ii A, B and C are as yet unknown functions of r, cr^ and cr^. We 
assume that ^o? he in the allowed range of the ^ coordinate - roughly 
speaking this is true when the width of the domain wall is small in comparison 
with radia of curvature of the core in the local rest-frame of the considered 
piece of the core. 

Notice that ±$o are exact solutions of Eq.(2). They are defined also 
outside of the region of validity of the co-moving coordinates. Therefore, 
formula (18) actually defines the field $ in the whole space-time. 

The Ansatz (18) implies that in fact we introduce boundaries of the do- 
main wall. The outer boundary in the co-moving reference frame is 
defined by the formula 

X(+)(r, a\ a^) = X(r, a\ a^) + ^o(r, a\ a^) n{r, a\ a^), 

while for the inner one (-^(-)) 

X(_)(r, a\ a') = X(r, a\ a') - ^(r, a\ a') n(r, a\ a'). 

Here n is the spatial part of the four- vector (n^). 

Inserting the cubic polynomial into Eq.(17) and equating to zero terms 
proportional to the zeroth and first powers of ^ we obtain the following 
recurrence relations: 

B = AKl (19) 
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C ^ -n''^^A+ {K^K^^-2X^l)A + K;^B, (20) 

where 

= -^i^(V^9"'4i:) (21) 

is the three-dimensional d'Alembertian on the world-volume E of the core. 
Of course, the cubic polynomial (18) does not obey Eq.(17) exactly. The 
leftover terms in Eq.(17) are of the order C,'^ and higher. Wc assume that 
these terms are not important. We will discuss the problem of accuracy of 
our approximation in the Appendix. 

We also require that the field $ is continuous everywhere, in particular 
at the boundaries, i.e. for C — Co^C = Then, by a standard reasoning, 
we deduce from Eq.(17) that also 9^$ is continuous at the boundaries. The 
second derivative 9|$ is not continuous at the boundary, in general. The 
conditions of continuity of $ and 9^$ at ^ = ^ and ^ = — give 

Co = 6, (22) 



A^-^, B = 0, C^-^. (23) 



The condition B — together with relation (19) gives the equation of 
motion for the core S: 

K = 0. (24) 
It can be shown that this equation is equivalent to Nambu-Goto equation 

□ (3)j^/. = (25) 

for a relativistic membrane. Thus, in the approximations we have made, the 
core can be regarded as the Nambu-Goto type relativistic membrane. 

Relations (20), (23) give also the following 2-1-1 dimensional, nonlinear 
wave equation 

+ (i^ - K^KDA - = (26) 

for 74(t, (7^, (T^) = A is proportional to the inverse width of the domain 
wall (2^) measured in the natural length unit /q, and it can be regarded 
as a scalar field defined on the core of the domain wall. Let us note here 
that 2(^0 is the width in the co-moving coordinates. Its transformation to the 
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lab-frame is not trivial: it includes Lorentz contraction and other changes. 



This transformation has been discussed in detail in the paper [jT^ in the case 
of cylindrical domain wall. 

Eq.(26) can be written in terms of dimensionless variables: introducing 
dimensionless f = t/Iq and a^'^ = a^'"^ /Iq , = d/df^ etc., we have 

U^'^A + - llKlKt)A - ^A" = 0. (27) 

Because g"''' as well as IqK^K^ are dimensionless, all coefficients and variables 
in Eq.(27) are dimensionless. 

Certain solutions of equation (27) for A can be found in an approximation 
scheme which is a generalization of a perturbative expansion proposed in 



1^ in the case of cylindrical and spherical domain walls. In that paper 
one can find a detailed motivation for the subsequent steps; in the present 
paper we will only briefiy describe the scheme. Our scheme probably gives 
only a certain class of solutions - there might be other solutions which can 
not be obtained in this way. Roughly speaking, the idea is to restrict our 
considerations to the cases such that the dimensionless extrinsic curvature 
l^Kj^K^ is close to zero, and to expand A in non-negative powers of it. In 
the zeroth order approximation Eq.(27) has the constant solution 

= 1. (28) 

There are also other solutions in this order, e.g., the ones having form of 
waves propagating along the core. Solutions of this latter type are time- 
dependent and have a characteristic frequency of oscillations > 1, where 1 is 
the (dimensionless) mass of the A field. Notice that this mass coincides with 
the mass of the $ field in the unit Iq^. We assume also that IqK^KI^ and its 
derivatives are smooth functions of f,cr^,o"^, and such that the derivatives 
di){llKl^K^) are much smaller than the function l^K^K^ itself. In this case 
the term IqK^K^A present in Eq.(27) generates a smooth, non-oscillating 
component N in A. We introduce the following perturbative Ansatz for A: 

A = l + n + N, (29) 

where Q denotes the oscillating component - the amplitude of these oscilla- 
tions is assumed to be of the order l^K^K^ . Equation (27) can be split into 
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two equations: 




The term O^^^N is by assumption regarded as small in comparison with 
the zeroth order contribution to N. Therefore, eq.(31) implies the following 
result for the lowest order contribution A^*^^) to N: 



Prom Eq.(30) we obtain the equation for the zeroth order oscillating con- 
tribution fi(^) to Q: 



Detailed form of this equation depends on the metric gab- Because the fre- 
quency of oscillations implied by Eq.(33) coincides with the mass of the scalar 
field (measured in the units Iq^), the presence of these oscillations is quite 
natural. The domain wall with the oscillating width can be regarded as a 
perturbation of a proper domain wall, which by definition does not have the 
oscillating component. 

3 Dynamics of the core in terms of local cur- 
vatures 

Equation (24), which governs evolution of the core of the domain wall, has 
rather abstract form. In order to bring its contents to surface, we will write 
it in appropriately chosen local coordinates r, a^, on the world- volume E 
of the core S. We shall call them the physical coordinates. 

As the T coordinate we take the one present in formula (4), without 
any further specification, while the coordinates r, a^, o"^ are now defined in a 
vicinity of an arbitrarily chosen point Xq{t) of the core 5 in a special manner 
explained below. For X = (X^,X^,X^) in formula (4) we write 



(32) 



(33) 



X{t, a\ a^) = Xo(r) + ei(r)ai + e2(r)a2 




(34) 
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where x denotes the vector product, and Ci, 62 are unit vectors specified 
below. For our purposes, the terms 0^^\a) (third order terms in cr\ cr^) do 
not have to be written exphcitely. It is clear that — — a'^ corresponds 
to the point Xq{t), and that the vectors X^„i,X^„2 tangent to S at Xo{t) are 
equal to ei(T), 62 (r), correspondingly. From conditions (6) we obtain that at 

Mr) 

— * 

rh 
n — —1= _ =, n — —1= =, (35) 



-{X^mf \jl-{mXof 
where 

m = ei X 6*2 (36) 

is the unit vector normal to 5* at Xq{t) and Xq = cLXq/cLt. The extrinsic 
curvature coefficients at Xq{t) arc equal to 

Xqu, c,m 



l-(mXo)2 yi-(mXo)2 

as it follows from the definition (7). The fact that K12 vanishes means that 
the tangent vectors ci, 6*2 have been chosen in directions tangent to the circles 
of main curvatures of the surface S at the point Xq{t). The parametrisation 
(34) is the most natural local parametrisation of the core from the viewpoint 
of geometry of surfaces in the 3-dimensional space. 

As for the choice of Xq{t) at different times r, there is a large freedom 
due to the reparametrisation invariance of the equation (24). In physical 
terms, this invariance means that the core S is a. "structureless" surface in 
the sense that translations along S are not observable. Natural choice is that 
Xq{t) is a smooth function of r. Apart from this, ^0(7") for different r's can 

be chosen almost arbitrarily. Even the condition X < 1 does not have to 
be imposed, because the core of the domain wall is merely a mathematical 
construct; in fact there is a numerical evidence that in a related case of 
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vortices the core can move with superluminal velocity, see e.g. |T^. We will 
use this large freedom to choose Xo{t) for different r in such a manner that 

lo(r)ei(r) =;eo(r)e2(r) = 0, (38) 

i.e. the motion of the point Xq{t) is always in the direction perpendicular 
to S. Taking (38) into account, we can write 

Xo{t) = v{T)m. (39) 
Then, simple calculations show that Eq.(24) reduces to the following relation 

*W + (40) 

Thus, locally, the acceleration i) of the piece of the core is determined by the 
main radia of curvature, up to the Lorentz factor. 

In the particular cases of a sphere (i?i = R2 = R, v = —R), and of 
a cylinder (i?i = cxo, R2 = R, v = —R), the usual spherical or cylindri- 
cal coordinates, respectively, have the properties required by formula (34) 
and conditions (38), and (40) coincides with equations considered in [[1^] . 
In a more general case, the explicit constructing of the special coordinates 
(cr^,(T^), such that formula (34) and condition (38) hold, might be difficult. 
Nevertheless, formula (40) is very helpful in qualitative analysis of motion of 
the core. 

If we drop the gauge condition (38) while still keeping (34), then instead 
of formula (40) we obtain 

V - (mel)(Xoel) - (mel)(Xoel) = 

[1 - (io)^ + (ioel)^]^ + [1 - (io)^ + (ioel)l^, (41) 

where v = Xorh is the perpendicular component of the velocity Xq. 

Formula (40) can be reinterpreted in terms of a surface tension. To show 
this we pass to the rest-frame of a small piece dS of the core. (Observe 
that formula (40) contains the laboratory frame quantities: the time = r, 

the transverse velocity v = Xofh, and the radia -Ri, i?2-) The core element 
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dS is parametrised by the parameters a^,a'^ introduced by formula (34): 
= = corresponds to the "center" Xo(r) of the element dS of the 
core. The radia i?i,i?2 are now the rest-frame curvature radia, denoted by 
i?]'^**, i?2^**. For simphcity, we restrict our considerations to the case when 

the piece dS of the core neither rotates nor is deformed at the chosen instant 
of time. Then = = and m = 0. Let us assume that on a piece dl of 
the boundary d{dS) of dS acts a force dF of the magnitude uj dl {dl = \dl\), 

— * 

tangent to the core and perpendicular to dl, directed to the outside of dS. 
Here a; is a constant. It is easy to see that 

dP = u N{a\a^) x dl, 

where 



— * 

is a unit vector perpendicular to the core at dl. In the last formula, the 
vectors X'^T*, X^^f are the vectors tangent to S at the point X{T,a^,a'^), 
calculated from the rest frame counterpart of formula (34), i.e. 

XT-^^ + ^t^^ ^T-^2 + ^,rn. (42) 

Notice that e^*^** = Cj - the vectors tangent to dS at Xq{t) - do not change 
under the boost in the perpendicular direction given by fh. The same is true 
for the parameters a^,a'^. Next step is to calculate the integral 



dF, 

d{dS) 

which gives the total force acting on dS. In the limit dS shrinking to the 
point Xq{t), the result for the integral is 

11 

frest = '^i^r^t + 

where \dS\ denotes the area of dS. (In fact, it is the dominating contribution 
only; we have neglected terms which vanish faster than \dS\.) In order to 



12 



find the corresponding force in the laboratory frame, we apply the boost with 
velocity v = \dXo/dt\ in the direction rh: fiab = {V^ — v'^^^^'^frest, i-e. 

Ill 

The laboratory frame curvature radia are related to the rest frame ones by 
the formula 

jyrest 

vi — f 

Notice that Ri are bigger than - this is due to the fact that in the 

laboratory frame the surface S is flattened in the direction of motion by 
Lorentz contraction. Another way to obtain this transformation law is to use 
the formulae (37) and the fact that Kii^K22 at = cr^ = are invariant 
under Lorentz boosts in the fh direction. Relativistic Newton equation of 
motion in the laboratory frame has the following form 

rf^Xo uj\dS\ I 1 

where s denotes the proper time. We have assumed that the rest-frame 
"mass" of the piece of the core comes entirely from the surface tension and 
is equal to cjIrfS*!. Equation (43) reduces to relation (40): one should relate 
the s variable to the laboratory time r {ds = Vl — v'^dr), to substitute 
dXo/ds — V fhjsJX — v^, and to use the assumption m = 0. If the piece 
dS rotates or is subjected to a deformation, the above presented reasoning 
should be generalized by introducing an appropriate stress tensor for the 
core. We will not dwell on this. 

At the end of the previous Section we have obtained the first correction to 
the inverse dimensionless width of the domain wall in the co-moving frame. In 
the absence of the oscillatory component it is equal to N'^^^ given by formula 
(32). The quantity KIK\ present in that formula can easily be calculated in 
the physical coordinates (cr^, o"^) defined by formula (34). We do not assume 
here the gauge (38). We obtain the following result: at the point Xq{j) 

^ 1 + /o' KlKl (44) 
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where 



(1 



,2^2 



mi 



^1 ^2 



Rf i?2 R1R2 



(45) 



Ri 



+ (m2 + 



Here we have introduced a short notation for components of the velocity and 
of the vector m: 



me,- 



-me, 



1,2; 



XqITI. 



(46) 



When the conditions (38) are satisfied, in formula (45) one should put Vi = 0. 

Let us present some implications of the formula (45). First, if each piece 
of the core is at rest at certain instant of time, then (at this moment) Vi = 




mj, and therefore formulae (44), (45) give 




(47) 



Thus, bending the domain wall is associated with making it thicker (if we do 
not excite the oscillatory component). 

If a piece of the core has a nonzero velocity at certain instant of time, then 
the terms in the square bracket present on the r.h.s. of formula (45) may 
compensate the effect of non-zero curvature because of the minus sign. To 
check this possibility, we have considered the following class of exact solutions 
of Eq.(24): 



x^, t) 



X 



(48) 



where li,l2 are constants obeying the condition II + I2 = 1, and / can be 
any smooth function. This plane wave type solution gives a membrane "lev- 
itating" over the {x^,x'^) plane on the altitude x^ = f{lix^ + — t), with 
an infinite plane wave of constant shape propagating as a whole along the 
membrane in the direction {hjh) with the velocity of light. The fact that 
such solutions exist can be deduced from a particular domain wall solution 



found in [|20[. Straightforward computations give for the solution (48) 



^ - ^ 

R2 ' Ri 



(1 + /'2)3/2 



/2' 



(49) 
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It follows from formula (45) that in this case N^^'^ = 0. Thus, in the leading 
order the presence of the plane wave travelling along the domain wall does 
not influence the width of the domain wall. 



4 The energy of the domain wall 

The total energy and energy density are very important physical character- 
istics of a solution of field equations, and therefore we would hke to calculate 
them for our domain wall solutions. The form of expression for the total 
energy depends on whether one integrates over the hypersurface of constant 
r, or over the hypersurface of constant laboratory time The construction 
of the conserved energy-momentum in the former case has been presented in 
Section V of the paper |l7]. In the present paper we consider the standard 



energy, i.e. the one obtained by integration over the hyperplane of constant 

Because the total energy of the field is constant during the motion of the 
domain wall, we can calculate it at an arbitrarily chosen laboratory frame 
time Of course, in order to really have the energy constant in time one 
should use the exact solution of the field equation (2). It will not come 
out exactly constant if we calculate it for an approximate, time dependent 
solution, in particular for the one given by formulae (18). 

The energy-momentum tensor in our model has the following components 

T^^ = d^<l>d''<^ + r]^"'L, (51) 
with L given by formula (1). The lab-frame energy E is given by the integral 

E = f d^xT"^^. (52) 



In this formula the field $, which is a scalar with respect to coordinate 
transformations, can be regarded as a function of the (r, cr^, cr^, ^) variables. 
With the help of a formula for differentiation of composite functions T^^ can 
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be written in the following form 



yo'^ = [(M-i)°«$,„ + n'^^,^]^ (53) 

where a — 0, 1,2; vP is the ji — Q component of the four- vector (n'*); the 
potential 

Vm^\{^'-^1)\ (54) 

and the matrix [M"-^] (the inverse of [Mab]) is given by formula (16). 

Also the volume element (Px is expressed by the co-moving coordinates, 



The integration over the ^ coordinate is effectively restricted to the interval 
(Co) ~Ci)) because in the vacuum = 0. Observe that on the hyperplane of 
constant the r variable in general becomes a function of ^ ; this is because 

T -|- ^n°(T, (7*) — x^ — constant. (56) 

Thus, constant x^ does not correspond to constant r, except for the core 
where ^ = and x^ — t. This complicates very much calculation of the 
integral over ^ because in general we do not know the explicit form of the r 
dependence - for this one would have to know explicit solutions of Eq.(24). 

There is a particular case in which the calculation of the energy is rela- 
tively simple: when each piece of the core of the domain wall is at rest at 
the initial time Mathematically, this means that the first derivative of X 
with respect to x° vanishes. Then, one can show that n° = 0, hence x^ — t 
for all ^. 

In the following part of this Section we will use the physical coordinates 
defined by formula (34) and the gauge condition (38). In these coordinates 
da^da"^ — \dS\, the area of the infinitesimal element dS of the core. For the 
core at rest v — Q and Cj = 0. This implies that 

d'x^{l-^){l-^)di\dSl (57) 
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and 



Here we have also used the relation (40). Therefore, formulae (57), (58) are 
valid for the Nambu-Goto domain walls only. 

The next step is to integrate T^^ over ^ in the interval {^o,—^i)- The 
integrand is a rational function of ^. There is no danger of a vanishing 
denominator, because we have assumed that the width of the domain wall 
(2^o) is smaller than the curvature radia Ri, otherwise the region of validity 
of the CO- moving coordinates would be too narrow to cover the whole width 
of the domain wall. 

For $ we take the approximate solution (18) with (22), (23) taken into 
account: for — < C ^ Co 

$(r,<7^<7^0 = ?<fo(f -^C)- (59) 



2 "^eo 3e, 



We may write 



_ $0 9^0 . . X . _ ^0 9^0 . . . .„„^ 

- ^/i(^), V{^) ^ ^$^/|(.) = ||/|(.), 

where z = ^/^o has values in the interval [1,-1], and the dimensionless func- 
tions fo{z), fi{z), f2{z) are given by the following formulae 

/o = —zil - z\ h - ^(1 - z% h - - \z^f - 1. (61) 

The energy E = J d^x is given by formula 



E = \^l i \dS\ I dz 



O / , 

^ J the core 



\ds\ y\ 
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+ 




(62) 



where, by the assumption, $,o/Ri ^ 1- Let us recall that this formula is valid 
for domain walls with the core at the instant rest and obeying the Nambu- 
Goto equation (24). Up to this point we have not used the approximate 
solutions for the half- width following from the Ansatz (29). 

Let us point out two particular consequences of formula (62). First, we 
see that all nonuniformities of the width of the domain wall increase the 
energy density because nonvanishing derivatives d^o/dcr'^ always give positive 
contribution to the integrand. Also r-dependence of increases the energy. 

Second, there is a curvature-dependent energy associated with the width 
of the domain wall which is present even in the case of constant width. It is 
given by the second term on the r.h.s. of formula (62) ( the one with /i, /2). 

For the solution (59) the dependence on z is explicit, and it is not difficult 
to perform the integration over z. Because the resulting formula is lengthy 
we will not present it here in the general case. Let us calculate the energy 
in the particular case, when the oscillatory component is absent, the half- 
width ,^0 is constant, and each piece of the core is at instant rest. Then, 
is related to the curvatures by formulae (44), (45), with v = Vi = mi = 0. 
Because the core is at the instant rest, the curvature radia have vanishing 
derivatives with respect to r, and therefore d^o/dr — 0. Formula (44) for 
^q/21q is approximate, so we shall calculate E with the same accuracy, i.e. to 
the second order in lo/Ri. We obtain the following result 





(63) 




(64) 
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Immediate consequence of formulae (63), (64) is that domain walls with 
nonzero extrinsic curvature can have smaller energy. This is because c\ — C2 
has come out positive. For instance, a straight infinite cylindrical domain wall 
with the core radius R (in this case I/-R1 = 0,1/-R2 = 1/-R) has smaller energy 
per unit area than the planar domain wall (for which I/-R1 = 0, 1/-R2 = 0). 
Thus, the domain wall prefers to have wrinkles. Let us recall that this has 
been found with the help of the approximate solution. It remains to be 
checked whether c\ — C2 is positive for the corresponding exact solution of 
Eq.(2). For this reason we do not claim that this is a proven result - it is 
just an indication, to be checked in another investigation. 

Another interesting problem, namely finding local minima of the energy 
e.g. for fixed area 15*1 of a compact core with a given genus, we also leave 
for a future investigation. 



5 The axially symmetrical toroidal domain 
wall 

In this Section we shall apply the presented formalism to a toroidal domain 
wall. Our motivation for doing this is that such domain walls are next to pla- 
nar, cylindrical or spherical ones with respect to complexity of their geometry, 
and to our best knowledge they have not been investigated as yet. Cylindri- 
cal or spherical domain walls have been considered in, e.g., pl| , p , |ll3| , [p!^ 



We shall consider the simplest case of a toroidal domain wall, characterised 
by the axial symmetry with respect to rotations around the x^-axis. 

We shall parametrise the toroidal core of our domain wall (at each instant 
of time) by two angles, G [0, 27r] and d G [0, 27r]: 

/ [R{T)^r{r,e)cose) cos<p \ 
X(r,^,0)= [rIt) +rlT,e)cose) sincP , (65) 
\ r(r, 9) sinO J 

where r is the laboratory frame time as introduced by formula (4), (j) is the 
azimuthal angle in the {x^,x^) plane, and 6 is the angle parametrising cross 
sections C of the torus with the half-planes of constant </>. C does not have to 
be a circle. Because of the axial symmetry, the cross sections C are identical 
for all angles (p. The radius R{t) gives the distance from the x^-axis to a 



19 



"central" circle of the torus. Only if the cross sections C are circular there 
is a natural choice for R{t): the distance from the x^-axis to the centers of 
the circles C. In general, there is a freedom in the choice of the "central" 
circle; this circle is merely a mathematical construct - the physical object is 
the domain wall. The choice of R{t) has influence on the form of the r{r, 9) 
function. In the following we choose 

R{t) — constant = Rq. (66) 

Thus, we have to calculate only one function r(T, 6). For correctness of the 
parametrisation (65) it is required that r{T,9) > and Rq > r{T,9). It 
might happen that in order to follow evolution of the torus for a prolongated 
interval of time it is necessary to introduce several parametrisation patches 
given by (65) with different RqS. This is the case when after some time 
the "central" circle with the fixed radius Rq turns out to lie outside of the 
displaced torus. 

Inserting the Ansatz (65) into Eq.(24) we obtain after straightforward 
computations the following nonlinear equation for r(T, 0): 

r{r^ + r'^)r - 2rr'{f'r - rr') + (1 - r^)(r^ + 2r'^ - rr") 

where r = 9r(r, 9)/dT, r' = dr{r, 9)/ 89. 

In general, equation (67) is not identical with formula (40) - they are 

— * 

written in different coordinates. The point is that X given by formula (65), 
when expanded in 9 — 9Q,(f) — (f)Q in a vicinity of a point (^^O) 0o) , does not have 
the form (34), and also conditions (38) are not satisfied, in general. In other 
words, reparametrisation gauge fixing implied by the formulae (65), (66) is 
different from the gauge fixing implied by the formula (34) and condition 
(38). Equation (67) can of course be transformed to the form (40), but this 
transformation involves nontrivial changes of coordinates. 

For simplicity, in the following part of this Section we shall consider only 
tori which at certain initial instant tq have circular cross sections. Such tori, 
in addition to the axial symmetry have also a reflection symmetry — > —x^. 
We also assume that for r = tq each point of the torus has zero velocity. 
This case is sufficient in order to get a feeling about evolution of the toroidal 
domain walls. 
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One may ask whether the initially circular cross section will stay circular 
at later times. To investigate this, we recalculate the basic equation (24) 
for the Ansatz (65) without the assumption (66) that R is constant in time. 
If the cross section C stays circular, there should exist such a choice of the 
function R{t) that the solution r is constant in 6. We have found that such 
a choice is not possible. Therefore, the initially circular cross section will 
always be deformed during evolution of the torus. 

Qualitative picture of the motion governed by Eq.(67) can be obtained 
by considering motion of the outermost and innermost circles on the core, 
^ = and 6 = 71, respectively. Because of the symmetries of the torus, r' — 
for 9 = 0, TT, for all times. Therefore, for = Eq.(67) reduces to 

r ,1 r'\ 1 



1 — f2 r Rq + r 

Similarly, for 6* = tt we have 



(68) 



1 -2 V • (69) 

1 — r Kq — r 

One can check that equations (68), (69) coincide with formula (40). The 
main radia of curvature of the torus at the point {9, 0) are given by the 
formulae 

1 1 r'sin9 + rcos9 ^^^^ 



Ri yV^ + r^ Ro + rcos9 



(71) 



i?2 (r2 + r'2)3/2 ■ 

Because r' = for ^ = 0, tt, the curvature radia are, respectively, Ri = 
±(i?o ± r), i?2 ^ = 1/r — r"/r^ (we choose the convention that the normal 

m to the torus is directed inwards). Also X does not have the tangent 
components for ^ = 0, vr, and f = — r in the both cases. Now it is easy to 
check that Eqs.(68), (69) follow directly from formula (40). 

Equations (68), (69) imply that there are two classes of tori, differing by 
their motion. To see this, let us consider Eqs.(68),(69) at the initial time Tq. 
Then C is a circle. Let us denote its radius by Tq. It is clear that Tq = Tq = 0. 
Also, by the assumption, r = at the initial time. Therefore, on the r.h.s. 
of Eq.(68) we have the force 

1 1 
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which is always negative. It follows that the outermost circle "6 = 0,0 — 
variable" will always move towards the x^-axis. As for Eq.(69), the force is 
equal to 

1 1 

+ ^ ' 

and it is negative only for ro < Ro/2. For ro = -Rq it vanishes, while for 
ro > Rq/2 it is positive. Therefore, in this last case the circle 9 = tt will 
start to move away from the a;^-axis. To summarize, when the initial circular 
cross section C of the torus is small enough, the torus will shrink towards 
a circle in the (x^, a;^)-plane with the center on the x^-axis. On the other 
hand, if the initial circular cross section is sufficiently large, we expect that 
the torus will start to shrink towards the x^-axis. As for the intermediate 
case of ro = -Ro/2, initially the acceleration of the innermost circle = 
TT, — variable" vanishes. At later times, as the cross section C of the torus 
shrinks r" becomes negative (because r has now a local maximum at ^ = tt), 
and on the r.h.s. of Eq.(69) the negative term r"/rQ will appear. Therefore, 
we expect that in the case of initial data such that ro = Ro/2, vq = 0, the 
circle 9 — n will start to move away from the axis. Numerical solutions of 
equation (67) confirm this expectations, see Figs. la-=-lc. 

Till now we have discussed evolution of the core of the toroidal domain 
wall. From formulae (44), (45) we obtain approximate value of the width of 
the toroidal domain wall in the absence of the oscillatory component. This 
is especially simple for 9 — 0,Tr, because there m — and Vi — V2 — 0. We 
find that 



2lo 1-v 



\, r'\„ 1 1 



r^ r (i?o ± r{Ro ± r) 



(72) 



where the ± signs are for 9 = 0,7r, correspondingly. From these formulae one 
can see that for r/Ro small enough (a slim ring far away from the x'^-axis) 
the outer side {9 — 0) of the toroidal domain wall is thicker than the inner 
side {9 — n). On the other hand, if r/i?o becomes closer to 1 (a fat ring close 
to the x^-axis) then, vice versa, the inner side is thicker than the outer one 
(when the oscillatory component Q is absent). 

As for evolution of the symmetric toroidal domain wall at later times, 
when the rest-frame curvature radia become too small for the Ansatz (18) 
to be a reasonable approximation, one can make some guesses based on 
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results of numerical computations, We expect that the slim torus will 
completely disperse into radiation after the first collapse to a circle in the 
(x^, x^)-plane. The fat torus will first change its topology from toroidal to 
spherical one - this will happen when the innermost circle of the torus is 
shrunk to the origin = x"^ = = 0. In the process some energy will be 
lost into radiation. Next, the domain wall will evolve as a deformed sphere. 
Eventually it will collapse to the origin, and it will disperse into radiation. 
This scenario might be different if we pass to another model. As pointed 



out in ||T^, sine-Gordon domain walls can reemerge after the collapse. In 
this case the domain walls can pass through each other with only a partial 
loss of energy into the radiation, and therefore one can expect that the torus 
collapses and reemerges several times before it follows the previous scenario in 
a final collapse. It has also been noticed in the paper [|13| that cylindrical sine- 
Gordon domain walls do not bounce in contradistinction to the spherical ones. 
This suggests that whether the sine-Gordon torus bounces or disappears after 
the first collapse might also depend on the big radius of the torus. At the 
initial time it is equal to Rq, and the curvature 1/Ri depends on it, see 
formula (70). For large Rq the torus locally is like a straight cylinder, while 
for small Rq curvature is more pronounced and such torus might behave more 
like a spherical domain wall. Probably the only way to verify these scenarios 
is to perform numerical calculations. 



6 Remarks 

(a) Calculations presented in Section 2 of this paper show that the method 
proposed in fl^ for cylindrical and spherical domain walls can be applied 



also in the more general case. The resulting equation (24) for the core of the 
domain wall is equivalent to the Nambu-Goto equation (25). Our method 
seems to be quite universal. Change in the basic equation (2) would result 
only in different form of the reccurence relations (19), (20). For example, if 
we change the quartic potential V (formula (54)) to the sine-Gordon type 
then only numerical coefficients in (19), (20) are different. The polynomial 
approximation can also be applied to non-relativistic domain walls in con- 
densed matter physics. 

We have not obtained any terms with higher derivatives in the equation 
for the core. Such terms are usually related to so called extrinsic curvature 
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corrections in the effective action for tlie core, ||^. Nevertlieless, we would not 
conclude that the pure Nambu-Goto equation for the core is the exact, final 
answer for all domain walls with small curvature, for two reasons. First, 
our solution is approximate one, and a better approximation could reveal 
corrections to the Nambu-Goto equation. Second, the basic field equation 
(2) possesses infinitely many solutions in the topological class of the single 
domain wall. Our cubic polynomial Ansatz probably picks only a subset of 
these solutions. There might be other solutions for which the core would not 
obey the pure Nambu-Goto equation (25). 

Actually, it is easy to point out possible ways to improve our approxima- 
tion and to find more general domain wall solutions. Natural step to find a 
wider class of solutions within the cubic polynomial approximation is to di- 
vide the interval (^o; (see formula (18)) into two subintervals, and to use 
two independent cubic polynomials in each of them. Next, one should match 
smoothly the two polynomials with each other and with the vacuum fields 
to obtain continuous $ and 9^$. In order to improve our approximation 
one could use polynomials of higher order, as suggested by considerations 
presented in the Appendix. Work along these lines is in progress. 

(b) Formula (40) from Section 3, giving the transverse acceleration in 
terms of the local curvature radia, is very useful in providing qualitative un- 
derstanding of time evolution of the core. For example, one could foresee 
from it that a bump on otherwise almost flat core will first split into two 
bumps travelling in the opposite directions, next more smaller bumps will 
appear and finally the bump will disperse into small riples. Numerical com- 
putations we have carried out in the case of an axially symmetric bump on 
the toroidal core described by equation (67) confirm this expectation. 

(c) Let us end this Section by mentioning several possibilities to extend 
our work. The first one is to improve our description of dynamics of the 
domain walls by including possible radiation, creation of pairs " domain wall 
+ anti-domain wall", and interactions between close parts of the domain 
wall. 

Second, one would like to have an analytical (approximate) description 
of domain walls also in the case when the curvature radia of the core are 
comparable with the width of the domain wall. Our results are for the case 
when the radia are much larger than the width. 

Third, one could try to apply the polynomial approximation to curved, 
non-static vortices in relativistic field theories as well as in condensed matter 
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physics. We have seen that this approximation is capable of yielding rather 
detailed information about dynamics of the width of the domain wall. It 
would be very interesting to have such an information also for vortices. 



THE APPENDIX. The accuracy of the polynomial approxima- 
tion 

Here we would like to comment on accuracy of the approximation used 
in this paper. It is instructive to apply the above presented scheme to the 
planar domain wall given by the exact solution (3). Then, 

Kab = 0, (gab) = diag {-1,1,1), ^ = x^, 

and in the absence of the oscillating component A = 1, i.e. = 2/o- Then 

$(x^) = $«(x^) = |$of (l-^(fr). (73) 

The boundaries of this domain wall are at = ±2/o. Comparing ^^^^ with 
the solution (3), we see that it is not equal to the first two terms of Taylor 
expansion of the solution (3). Thus, we do not recover the exact solution (3) 
term by term in the Taylor expansion. On the other hand, the approximate 
solution has the right global characteristics of the domain wall like energy 
per unit area or the boundary conditions. The energy per unit area in the 
case of solution (3) is 

6 to 



It is minimal in the topological class of single domain wall, ^2[. The energy 
of the solution ^^^^ is a little bit higher 

^(1) ~ 0.79^$^ 



0- 



For a comparison of energy densities, see Figs.2a,b. 

Improving the approximation (73) consists of including terms of higher 
order in x^. We know that the exact solution is an odd function of x^, so the 
first term to be included is proportional to the fifth power of x^. However, 
it turns out that a fifth order polynomial in x^ can not simultaneously obey 
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recurrence relations obtained from Eq.(2) and the continuity conditions for 
$ and (9a;3$. The reason is that the recurrence relations imply that the fifth 
power of comes with a positive coefficient; then it is easy to see that there 
can be a problem with smooth matching of the polynomial with the constant 
vacuum solutions. 

Going to the seventh order polynomial does yield a solution. Denoting 
z = x^/lo, the result we have obtained can be written in the form 

^ $0 (0.5750Z - 0.0479/ + 0.0059/ - 0.0004/). 

smoothly matches the vacuum solutions for ^ ±2.593Zo, while for the 
solution the matching takes place at x^ — ±2lo. Energy per unit area 
for the solution ^^^^ is equal to 

fti 0.706^$^, 

which higher than Eq by only 6%. Also the energy densities do not differ 
much from their exact values, see Figs.3a,b. 

To summarize, the direct comparison presented above in the case of the 
planar domain wall shows that the polynomial approximation works rather 
well. Therefore, we expect that also for slightly curved domain walls, i.e. such 
that their rest-frame curvature radia are large in comparison with their width, 
the cubic polynomial approximation is a reasonable one, and that in order to 
improve it one could just use a higher order polynomial. Calculations with 
seventh order polynomials do not have to be cumbersome, because Taylor 
expansions and other necessary operations can be carried out by a computer. 
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Figure captions 



Fig. la. Numerical solution of equation (67). The initial data are r{0,9) — 
0.4, r(0, 9) = 0. The contours are the cross sections C of the torus with, e.g., 
{x^,x'^) plane. They are determined by the function r(t,6) (with 6 changing 
from to 27r). The cross sections are shown for the following values of 
t = t/Rq: 0.00, 0.15, 0.25, 0.35, 0.45, 0.55. The horizontal axis shows the 
distance from the axis, measured in the unit Rq. 

Fig. lb. The same as in Fig.la, but with the initial data r{0,6) = 0.5, 
r{0,e) = 0. The values of t are: 0.00, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65. 

Fig.lc. The same as in Fig.la. but with the initial data r{0,6) = 0.65, 
f{0,e) = 0. The values of t are: 0.00, 0.15, 0.25, 0.35, 0.45, 0.55, 0.65. 

Fig. 2a. The comparison of the gradient contributions l/2{dxS^y to the 
energy density of the planar domain wall. The solid line corresponds to the 
approximate solution ^^^^ , the dashed line corresponds to the exact solution 
(3). 

Fig. 2b. The comparison of the potential energies V{^) (formula (54)) for the 
planar domain wall. The solid line corresponds to the approximate solution 
the dashed line corresponds to the exact solution (3). 

Fig. 3a. The same as in Fig. 2a, but for the approximate solution 

Fig. 3b. The same as in Fig.2b, but for the approximate solution 
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